MEROMORPHIC SOLUTIONS TO A DIFFERENTIAL-DIFFERENCE 



EQUATION DESCRIBING 
CERTAIN SELF-SIMILAR POTENTIALS 

ALEXANDER TOVBIS 

Abstract. In this paper we prove the existence of meromorphic solutions to a nonlinear dif- 
ferential difference equation that describe certain self-similar potentials for the Schroedinger 
operator. 



1. Introduction 

Let L = —d 2 + u(x) be a Schrodinger operator that we factorize as L = A + A" + A, 
where A^ 1 = ^fd x + f(x) and x, A G C. Then the function f(x) satisfies the Riccati equation 
f 2 (x) — f'(x) + A = u(x). If L = A~A + + A denotes a new Schroedinger operator, obtained 
from L by permuting the operator factors A ± , then the potential u(x) of L is given by 
u{x) = f 2 (x) + f'(x) + A. The differential-difference equation (DDE) 

(1) [f(x) + f{x + a)]' + fix) - fix + a) — fj,, 



where a and fi are some complex constants, was derived in [Si], to describe potentials of 



the Schrodinger operator L satisfying the self-similarity constraint L = ULU 1 + /i, where 
U is the translation operator U fix) — f(x + a). In the case /i = equation ([J) is satisfied 
by(0,H) 

_ l V'ix-x Q )-V'ia) 
[) J[) ~ 2P(x-x„)-?(a)' 

where V denotes the Weierstrass elliptic function and xq G C is an arbitrary constant. The 
function (fj) is a meromorphic function with only simple poles. 
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The aim of the present paper is to prove existence of meromorphic solutions with simple 
poles to (P in the case /i 7^ 0. Without any loss of generality, we assume a to be a positive 
real number. Indeed, if a = re^, where (f> G R and r > 0, then the transformation x — ► e^x, 
f — > e~^f reduces ([[]) to an equation of the same type with the step r G R + . Our approach, 
which utilizes some ideas of |To|, consists of essentially three statements: 1) there exist two 
different formal power series solutions (in powers of x~^) to ([[]); 2) for any formal solution 
f(x) there exists an actual solution f(x), analytic (when \x\ is sufficiently large) in a sector S 
on the complex x-plane of opening greater than tv and having the asymptotic expansion f(x) 
in S; 3) any such solution can be meromorphically continued onto C so that f(x) may have 
only first order poles. These statements are proven in Sections [3] and ^] respectively. The 
author wants express his gratitude to V. Spiridonov for interesting discussions held at the 
NATO ASI "Special functions 2000" , Tempe, Arizona and for the following correspondence. 
The author also want to use this opportunity to thank the organizers of the NATO ASI 
"Special functions 2000" . 



2. Meromorphic continuation 

Let R G C be a simply-connected domain, bounded by piece-wise smooth curves 771 (0 
and 772(C)) where x = C + in is a complex number and 772(C) > 771(C) for all £ G (— 00, +00). 
The values r] 2 = +00 and 771 = —00 are allowed. We say that the domain R is a-wide on the 
interval / C R if r/ 2 (£) - 771(C) > a for all (el. 

Theorem 2.1. A solution f(x) to the equation (|I|) that is analytic in some domain R that is 
a-wide on R admits a meromorphic continuation on the whole complex plane and all possible 
singularities of f(x) are first order poles. 

Proof. Let Ri = R and let R n denote the domain bounded by the curves 771(C) and 772(C) + 
(n — l)a, where n = 2,3,.... We first prove that f(x) can be meromorphically continued 
to the domain Roo = UfR n , i.e., to the right of the domain R, by considering (|]) as the 
Riccati equation 

(3) f(x) = f(x) + h(x), 
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where h(x) — /j, — f(x — a)' — f 2 (x — a). The latter equation is equivalent to to the second 
order linear differential equation 

(4) u"{x) + h(x)u{x) = 
through the standard transformations 

(5) u{x) = exp / f(t)dt and f(x) - 

J X 



u'(x) 



u(x) 

where x G R is a point such that x — a G R. Note that, according to the function u(x) 
is analytic in R. 

Consider (f|) in the region Qi, where Q n = R n+ i\R n , n G Z + . According to the assumption 
of the theorem, the function h(x) is analytic in Qi, so that solutions of the linear equation 
(^) are analytic in Q\. Thus, we can analytically continue u(x) onto R2. If u(x) does not 
attain zero value in Qi, we get an analytic continuation of f(x) onto R2 by (|5|). However, 
if u(xq) = for some x G Q±, then f(x) has a first order pole in Xq, which is an isolated 
singularity of f(x). Thus, we obtained the required meromorphic continuation of f(x) onto 
i?2- This process can be continued to the domains R%, R4, etc. in the same fashion. However, 
now we have to consider a possibility that h(x) has a singularity at x + a. Then, u(x) and, 
correspondingly, f(x) may have singularities at x n = xq + na, where n G Z + . We need to 
show that these possible singularities of f(x) are first order poles only. 

Let 

h(x) = h + hi(x — x ) + h 2 (x — x ) 2 + . . . and 

(6) u(x) — (x — Xq) + U 2 (x — Xq) 2 + U 3 (x — Xq) 3 + . . . 

be the Taylor expansions of h(x) and u(x) near x — Xq in (§]). Comparing the like powers of 
x — xq in ([|), we obtain 

(7) u 2 = 0, u 3 = - — , Ui = -—, 

o 12 

so that the principle part of f(x) at x = xq is —{x — xq)^ 1 . 
Combining the expression for h(x) with (Q), we obtain 

(8) + a) = /x - - 2f 2 (x) = 11- h{x) - 2[(lnu(x))'] 2 . 
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Direct computations show that 

(9) [\nu(x)}' = — ^(x - x ) - ^(x - x ) 2 + 0(x - x ) 3 , 

x — x 3 4 

and 

(10) ~2([lnu(x)]') 2 = — 2 + ^ + h(x - x ) + 0(x - x ) 2 , 

[X — Xq) 6 

so that 

—2 h 

(11) h(x + a) = h(x + a) = - + (// + -^) + 0(x - x ) 2 . 

Considering now this equation near the point x 1 e Q 2 , we obtain 

(12) h(x) = ~ 2 + (/x + ^) + 0(x - x x ) 2 . 

The proof that all the points x n are either regular points or first order poles of f(x) follows 
by induction from the following two lemmas. 



Lemma 2.2. The differential equation 

(13) u"{x) + h n (x)u{x) = , 

where the coefficient h n (x) has the form 
(14) 

h n (x) = ~ -^—^ +A + A 2 [x - x n f + ■■■ + A 2N _ 2 (x - x n ) 2N ~ 2 + 0(x - x n ) 2N ~ l , 

yX X n ) 

posseses two linearly independent solutions 

(15) u(x) = (x- x n ) N [I + u 2 {x - x n ) 2 + ■■■ + u 2N {x - x n ) 2N + 0(x - x n ) 2N+l ] 

(16) v(x) = {x- x n )~ N+1 [I + v 2 {x - x n f + ■■■ + u 2N {x - x n ) 2N + 0(x - x n ) 2N+1 ] 

that are analytic in a neighborhood of x = x n . Here Ak,iik,vt G C and n, N e N with 
n + 1 > N > 2. 
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Proof. Frobenius multipliers of (|T3D at N and —N + l, so the first terms of u(x) 

and v(x) are (x — x n ) N and (x — x n )~ N+1 respectively. Suppose, 

oo 

(17) u{x) = (x - x n f + Y^u k (x- x n ) N+k . 

k=l 

Computing w fc , we see that the odd coefficients u± — u 3 — ■ • ■ — u 2 n~i = 0, so that u(x) is 
in the form fll5|) . Similar arguments work for the second solution v(x). □ 



Arguments of Lemma 2J2 show that solutions to (|13[) have no branching at x = x n since 



Frobenius multipliers are integer and there are no logarithms. Note that the general solution 
to ( O ) has the form 

(18) u(x) = (x- x n )~ N+1 [u + u 2 {x - x n f + ■■■ + u 2N _ 2 {x - x n ) 2N - 2 + 0(x - Xn) 2 ^ 1 ] 

where u ^ 0. The only nontrivial special solution to (|13|) is proportional to u(x) given by 
flT5|). In any case, the point x n is a simple pole of the solution f{x) = —\\nu(x)]' to (H). The 



following Lemma shows that the new coefficient 

(19) h n+ \(x) =ii— h n (x — a) — 2[(\nu(x — a))'] 2 , 



where u(x) is a solution to (0), has the form QliD with another N. Thus, the point x n+ \ is 
also a simple pole or a regular point of the solution f(x). (Note than x n+1 is a regular point 
of h n+ i if the new N + l.) 

Lemma 2.3. The new coefficient h n +i(x) has the form ( fT^ ) with the new N equal to N — 1 
if u(x) is given by (|I8"D or equal to N + 1 if u(x) is given by ([15]) . 



Proof. Consider, for example, the case when u(x) is given by (J15|) - Using expansion (|14D for 
h n (x), we obtain the first odd coeficient u 2 ^ + i in ( ]T5| ) is 

/ on N A 2N+ i 
(20) M27V+1 " _ 4iV(2Ar+l) • 

On the other hand, we see that hau(x) = A^ln(x — x n ) + u 2 (x — x n ) 2 + . . . , where the 
leading odd term of \nu(x) is u 2 ^ + i(x — x n ) 2N+1 . Then 

[\nu(x)}' = {x- x n y l [N + 2u 2 (x - x n ) 2 + . . .] 
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where the leading odd term in the square brackets is (2N + 1)m 2 at+i(x — x n ) 2N+1 . Finally, we 
get the leading odd term of —2([\nu(x)}') 2 as — 4N(2N + 1)u2n+i( x — %n) 2N ~ l ■ So, according 
to fli~9f) and (p0|), the leading odd term of h n+ i(x) at x n+ i is of the order (x — x n+ i) 2N+1 . The 



leading term of h n+1 (x) is ^~~^r = — (^z^~y2 ■ So, /i n+ i(a;) has the form ([TJ ) with the 



new iV equal to A" + 1. The case when u(x) is given by (|i"8|) can be considered in a similar 
way □ 

Consider, for example, the singular point X\. According to (|T2|), we have N = 2 ah this 
point. Solving the corresponding initial value problem for (|l"3"D, we get the solution u(x) of 
either type fll5D or ([181) . In any case, f(x) has a first order pole. The corresponding function 
h,2 (x) , according to Lemma has a second order pole at x = X2 with the principal part 
— ( x _^ 2 )a if u(x) is proportional to ([15|), and is regular at x = x<i if is given by (plf). 
We can continue these arguments to show that all points x n are either first order poles or 
regular points. Thus, we proved meromorphic continuation of f(x) on R^. 

To prove the meromorphic continuation to the left of the domain R let us note that the 
transformation x = — t — a, g(t) = f(—t) reduces the equation ([]]) to the equation of the 
same type 

(21) [g(t) + g{t + a)}' + g 2 {t) - g 2 (t + a) = - M , 

which have an analytic solution g(t) on t G —Roo- We can now use the previous arguments 
to continue g(t) to the right on the whole complex plane. □ 

Corollary 2.4. Let I cK be an interval and let f(x) be a solution to the equation ([!]) that 
is analytic in a domain R that is a-wide on I. Then f(x) admits a meromorphic continuation 
on the strip Qx G / and all possible singularities of f(x) are first order poles. 

3. Asymptotic solutions 
In this section we rewrite ([[]) as 

(22) nf(z)=fi + Af 2 (z) , 

where z = x + b, b = f and the operators □, A act on a function g(z) as 

(23) Og(z) = g(z + b) + g(z - b) and Ag(z) = g(z + b) - g(z - b) . 
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Some important for us properties of □, A are given by identities 

(24) 2D[fg] = □/•□# + Af ■ Ag and 2A[fg] = Af ■ Dg + Df ■ Ag , 

where /, g are given functions. 

In this section we construct two different formal power series solutions f ± {z) to (fl|) and 
prove existence of corresponding actual solutions (z). These solutions are analytic in the 
corresponding sectors S^, n = 0,1, 2, specified below, when \z\ is sufficiently large and admit 
asymptotic expansions 

(25) z^oo, zeS±. 



This fact together with Theorem ^Tj prove existence of non-trivial meromorphic solutions to 
©)■ 

Proposition 3.1. Equation (^) possesses two formal power series solutions 
(26) /±(*)=±AVi+2£ + ZX*-* , 

fc=2 

where A = ^ and i/ie coefficients y^ are defined uniquely. 



Proof. The substitution / = \zz + c + y(z) reduces the equation ( 22|) to 



(27) ^d* - * +Dy' = fi + (XAz^ + A#)(AOz* + 2c + Dy). 

Comparing leading coefficients and taking into account Az^Ozi = 2b, we get A = ±^/— ^. 
Taking into account (]24]) , the latter equation can be now rewritten as 



(28) 2AA[^£] = □y / - Ay 2 - 2cAy + A(-Dz"5 - 2cAz?) . 

To expand y in powers of z~^ we need the free term of ( ^8"D to be of order 0(z~ r ) or less. 
Thus, we obtain c = ^, so that the free term is of the order 0(z~i). Equation (28) can be 
now rewritten as 

(29) A [z -,y] — — (Oy - Ay - ^Ayj H - — . 
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It is clear that the expression in the left hand side is the dominant term of the latter 
equation and that the substitution y = Vk z ~^ satisfying (p9|), defines the coefficients 

y k uniquely. □ 

Let f(z) = X^fz + + • • • denote one of the formal solutions f^(z). Given A G C, we 
define angles 

(30) /? n = ^7r(l + 2n)-^argA, 

where n = 0,1,2 and argA G [0,2-71"). If (3 n + tt is not a multiple of 27r, we define the 
sector 5+ on the Riemann surface of z 1 ^ 3 by extending a small sector bisected by arg z = (3 n 
independently in both positive and negative directions until it either hits the negative real 
direction or the ray axgz = f3 n ± tt respectively. If j3 n is not a multiple of 27r, the sector S~ 
on the Riemann surface of z x l 3 is defined by extending a small sector bisected by arg z = (3 n 
independently in both positive and negative directions until it either hits the positive real 
direction or the ray argz = (3 n ± tt respectively. In the case (3 n = 2nk, k G Z or f3 n = 
(2k + 1)71, k G Z, the sectors S~ or respectively are considered to be empty. 

Theorem 3.2. If f(z) is a formal solution to the equation ([!]) and is a nonempty sector 
described above, then there exists and actual solution f(z) to ([!]) that is analytic in sufficiently 
remote part of any proper subsector S of S£ and 

(31) f(z) ~ f(z) , z^ oo, z G S. 

Proof. Our main idea is to reduce the considered DDE to an integro-differential equation 
(IDE) and to show that the latter equation can be solved by successive iterations in a proper 
sectorial neighborhood of infiniti. Substituting f(z) = X\/z + ^ + y(z), we obtain equation 
(p^) for y(z). The inverse Laplace transform C" 1 , applied to (p9|), yields (see, for example, 
[Ob]) 

1. 



— 2smh(pb)C 1 [z 2 y](p) = — ( — 2p cosh(pb)Y (p) + 2sinh(6p) 

2A V 

sinhte. , , 1 . 
32 +-— Hcoth&p--) , 

2UTvp bp 



Y* 2 (p) + ^Y(p) 
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where Y{p) = £ _1 [?/](p) and F{p)* 2 = J Q P F(p — q)F(q)dq. After a simple algebra, the latter 
equation becomes 

(33) CrVM-myU-^)-^, 

where = xcothx — 1. Note that is a meromorphic function with simple poles at 
the points ink, where k G Z\{0}, and that H(x) has not more than linear growth in any 
non- vertical direction argx = const. 

Separating the linear part bp of 5(6p) along the positive real axis, we can rewrite (^) as 

(34) c-^m - irw - - iy-Op) - ^ . 

where the function S(fep) — bp is bounded on [0, oo). Applying now the Laplace transform L 
to d§, we reduce @3) to the IDE 



(35) 3/(2) + 2Az*y(z) = h{z) * y{z) - y\z) - T^=r{z) , 

where £(z) = C\E{bp) — bp](z), r(z) = C[E(bp)p^ 3 ^ 2 }(z) and the convolution is defined by 

i pA+ioo 

(36) f{z)*g{z) = —\ f(s)g(z-s)ds 

with a sufficiently large A > 0. Considering (|35"D as a perturbed linear equation y'(z) + 
2Xz^y(z) = 0, we rewrite the former as 

(37) y(z) = e~l Xz3/2 [ e^ /2 W[y{t)]dt , 



where the nonlinear operator W[y] denotes the right hand side of (|35|) and the contour of 
integration j(z) is to be specified below. Equation ( |3"7D can be rewritten in the operator 
form as y — TW[y). 

Let us assume for a while that | < argA < and consider sector Sq. According to (|30|), 
this choice of A allows us to find a proper closed subsector S C Sq that contains the right 
half-plane. We want to solve the IDE fl35| ) in a sufficiently remote part of the sector S by 
successive approximations. In order to formulate the statement more precisely we need to 
introduce the following notations. 
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Let S be the image of the sector S under the transformation £ = z 3 ^ 2 . Let zq be a 
sufficiently remote point on the ray argz = /3q and let S^ , where Co = z^ 2 denote the 
parallel shift of S so that the vertex of £ is shifted to Co- For every C £ E^ we define a 
contour r(£) as a ray eminating from £ and such that: T(Q = {r : arg(r — C) — f A)} 
if |argC - §A)| < f; T(C) = {r : arg(r - C) = argC + f} if argC - |/?o < -f and; 
F(C) = { r : arg(r — C) — arg £ — | } if arg £ — |/? > f • Then the region S ZQ and the contour 
j(z) are the images of E^ and T(Q under the transformation z = £ 2//3 . 

Let y (z) = 0, y k (z) = I o fc = 1,2,... and % = y k - y h _ x . We will show 

that the solution to fl3?D is given by 

oo 

(38) y(z) = X>k, 

fc=l 

where the series converges absolutely and uniformly in S Zo for sufficiently large \zq\. This can 
be done by introducing the Banach space B of functions h(z), such that h is analytic on S ZQ 
and satisfy \h(z)\ < B\z~ 2 \ there with some constant B > 0. (Note that B depends on h.) 

The norm of h £ £> is sup 2g5 \z~ 2 h(z)\. According to Lemma 14.2 from |Wa|| , the integral 

_ i 

operator I : B — > i3 is a bounded linear operator, where ||X|| is proportional to |z 2 1- We 
start to study the nonlinear operator W[y] by considering the convolution 



£( z ) * y( z ) = £[(S(6p) - 6p)y(p)](z) = £[{£>p(coth&p - 1) - e - 2fep }y(p)](z) 
(39) + C[(e~ 2bp - l)Y(p)](z) = n(z) * y(z) + y(z + 2b) - y(z) , 



where, according to ||OB 



(40) fa) = C[{bp(cothbp - 1) - e" 2 *] = + |) - — . 

Here ^(x) = Jj^y denotes the logarithmic derivative of the Gamma-function. It is clear that 
fi(z) is a meromorphic function with double poles at z — —2kb, k £ Z + . The asymptotic 
expansion 

(41) *'(*) ~ I + ± + £ ^ 



fc=l 
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where B 2 k are Bernoulli numbers, follows from the Stirling formula (see ||GR]| ) . Combining 
the latter facts, we obtain the estimate 

(42) Hz)\< zeS, 
where the constant \l/o depends only on the sector S. 
Proposition 3.3. If y E B then fi(z) * y(z) E B and 

(43) ||^*y||<M||y|| , 
where M > does not depend on y and zq. 

Proof. Consider first the case z belongs to the right half-plane S = {z : 5fc > mz }, where 
the number m > is choosen so that S C S Zo . Note that m does not depend on \zq\. Setting 
A = mz in (|36l) , we obtain 

pA+ioo 

(44) n(z)*y(z) = — y(s)fx(z- s)ds 

*m J A -ioo 

Note that we are integrating over the boundary of S and thus $l(z — s) > 0. The fact that 
fi(z)*y(z) is analytic in S follows from the properties of /i(z) and y(z) immediately. Utilizing 
and the fact that y E B, we obtain 

— T /i (z) * y{z)\ < — / 

/" +0 ° ^ 

27 r y_ 00 (A 2 + 7 ? 2 )[( U -A) 2 + (^-7 ? ) 2 ] ' 

where s = A + in and 2 + 26 = u + zt> with r),u,v G M. 

Let / be the latter integral and p(r/) denote its integrand, which has simple poles at the 
points rj = %A and n = v + i(u — A) in the upper half-plane. Computing I via the residues 
of g(rj) in the upper half-plane, we obtain 

1 



(45) 



/ = 27ri Res g{rj) 



TV 



v + i(u — 2a) 



+ 



A{y — iu) (u — A) (v + iu) 



(46) = A^TW-{ 1 + ^A 
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Thus 

< 47 > m^ )tm -m 

since A can be taken greater than 2b. 

Consider now the half-plane that is obtained by rotating S on the angle <p, where is 
choosen in such a way that S<f, C S ZQ . Define now another convolution m*^y by (f44[) , where 
the contour of integration is the boundary of S^. Using the same arguments as above, we 
can obtain the estimate 

(48) bM\ Hz)HyiZ)l -W' 

where z G and continuously depends on (p. However, m(z) y(z) is an analytic 
continuation of m(z) * y(z), since the functions coincide on S$ D S. Taking M = maxM^, 
where M — |^ we complete the proof of the proposition. □ 



To complete the proof of the theorem we use the standard technique to show the conver- 
gence of iterations (0). Using properties of the Laplace transform and of the operator I, 
one can show that y± G B. Let K = \\yi\\, and let us prove by induction that \\Sy n \\ < 2 1 ~ n K 
if \z 1 is sufficiently large. Indeed, according to the estimate of the convolution, 



/ AK 

(49) \\6y n \\ < IIXH IM + 2 + — 

where the induction assumption 



/cr\\ II » > li - li II \ S Vn-l\\ , „|(fy„-i| 



was used to estimate the nonlinear term of W. It remains to choose z so that ||X|| + 2 + < 
\ for z G to complete the proof of the theorem for | < arg A < -j 1 the sector Sq. 

Let us now consider the general case < arg A < 2tx. The sector Sq is given by — 77+2 A < 
argz < 7T if arg A G [0,7r] and by — 7r < arg ,2 < 7r ~ 2 ^ rgA if arg A G [7r, 27r]. Note that the 
opening of the sector Sq is greater than n for any arg A and that Sq contains the right hlaf- 
plane if ~ < arg A < ^f. Let us now choose a proper closed subsector S C Sq of opening 
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greater than tt, and let argz = a be the bisector of S. Clearly, \a\ < |. Let 

/•e~' a oo 

(51) C a [Y]{z) = / e- z "Y{p)dp 

Jo 

define the Laplace transform along the ray argp = —a. The contour for the corresponding 
inverse Laplace transform as well as for the corresponding convolution * a is a straight line 
perpendicular to argz = a. Therefore, we can use our previous arguments to show the 
uniform and absolute convergence of iterations (|38|) in a properly constructed S Zo C S. In 
the same fashion the theorem can be proven for any nonempty sector S£, tiSZ. 

Recall that the function H(6p) has poles on the imaginary axis and, therefore, £ Q S is not 
defined for a = ±~. However, we can define C a E for a such that \a — 7r| < ~ and repeat the 
previous arguments for sectors S 1 ", n e Z. 

□ 

Corollary 3.4. Let , n = 0,1,2 denote a solution of (HI), analytic in some remote part 
of S„ as described in \3. 1 . Then, according to \2.1\ is a meromorphic function on C t/jai 
can nave on/y simple poles. 
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